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Competition within the immune system may degrade immune control of viral infections. We 
formalize the evolution that occurs in both HIV-1 and the immune system quasispecies. Inclusion 
of competition in the immune system leads to a novel balance between the immune response and 
HIV-1, in which the eventual outcome is HIV-1 escape rather than control. The analytical model 
reproduces the three stages of HIV-1 infection. We propose a vaccine regimen that may be able to 
reduce competition between T cells, potentially eliminating the third stage of HIV-1. 
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Our immune system is highly effective in suppressing 
most viral infections, due to the many different T cells 
that exist in the repertoire of one person. While many 
different T cells can recognize a virus, only those of high- 
est affinity respond in large numbers and participate in 
eliminating the virus. One limitation of our immune sys- 
tem stems from the competition among T cells of similar 
specificity for the virus. For a series of discrete infec- 
tions over time, competition in the immune system is 
associated with the phenomenon of deceptive imprinting 
0, or original antigenic sin, and has been characterized 
by a random energy model 0- Original antigenic sin is 
the tendency for memory immune cells produced in re- 
sponse to a first viral infection to suppress the creation of 
new immune cells in response to a second infection with 
a related strain. Moreover, while they are used, these 
memory immune cells may not be optimal for control of 
this second, different viral strain. Another form of com- 
petition in the immune system occurs when several viral 
strains simultaneously infect one person. In this case, 
the T cells compete to recognize the different strains, and 
recognition of all strains may not be uniformly effective. 
This immuno dominance of one strain over others means 
that the immune response to multiple infections is not a 
simple superposition of the responses to each individual 
infection 

For a chronic infection with a mutating virus, the ef- 
fects of competition are more subtle. Since the infection 
changes slowly over months or years, there may be no 
apparent original antigenic sin, as the immune system 
slowly updates itself in an attempt to clear the virus. If 
the virus evolves mainly through point mutation, the mu- 
tants will be only slightly different from their parents and 
can still be eliminated by the memory T cells. Moreover, 
whenever a sufficiently new viral strain emerges through 
the accumulation of point mutations that cannot imme- 
diately be cleared by the immune system, there is enough 
time for the immune system to create new T cells against 
the new strain since the infection is chronic. Therefore, 
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FIG. 1: The entire time course of HIV-1 infection. The circles 
are clinical data from Q. The solid line is generated by our 
model. 



competition might represent a relatively small weakness 
of the immune system against chronic infections. How- 
ever, when facing HIV-1, any trivial weakness may be- 
come the immune system's Achilles' heel. A typical HIV- 
1 infection is characterized by three phases. See Fig. ^ 
The viral load initially increases dramatically in the first 
viremia, followed by a sharp decline due to immunologi- 
cal control. Thereafter, a long period of clinical latency 
lasts for 6-10 years with small viral load. Finally, the 
HIV-1 viral load climbs up again in the final viremia, in 
the onset of AIDS. Most modeling works focus on the first 
viremia (e.g., [ElI) In 0- the second and third stages 
were modeled but without the first stage. The entire time 
course has been difficult to reproduce faithfully by a uni- 
fied mathematical model, probably due to the disparate 
dynamics of the three phases. To describe each phase, 
a different set of parameters has been given. The three 
stage dynamics is reproduced by computational models 
m 01 > but the mapping between parameters and phys- 
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iology remains to be fully established in those models 

In this letter, we develop an analytical, physical the- 
ory of the immune response that explains HIV-1 patho- 
genesis. A key concept in the model is a population of 
viruses with different sequences that arise due to viral 
evolution, termed the viral quasispecies. Since HIV-1 
and the immune system are well-matched competitors, 
the viral quasispecies and its dual, the lymphocyte qua- 
sispecies, should be equally well matched. These con- 
siderations lead to a theory with a large number of dy- 
namical equations, which is reduced to a low dimensional 
model by symmetry. The model fits well with the stan- 
dard HIV-1 data of the entire infection time course, and 
the parameters are all biologically reasonable. A novel 
equilibrium between the immune system and virus qua- 
sispecies is found, which, instead of fixing the viral load, 
favors viral escape. 

We let Vi (i = 0, 1, • • • I - 1) denote the HIV-1 quasis- 
pecies i. We let Vi denote the concentration, or load, of 
Vi. We let v — J2i=o Vi denote the total viral load. We 
let Xi denote the T cell, or effector T lymphocyte, that 
is specific for Vi- We let Xi denote the concentration of 
this T cell, and we let x = X^=o Xi denote the total T 
cell concentration. We let kji denote the affinity between 
the T cell Xj and the virus Vi. As a fairly generic as- 
sumption, we take fey = exp(— fidij), where = dji is 
the Hamming distance between Vi and Vj . Since roughly 
two amino acid substitutions still permit recognition, we 
expect (5 to be on the order of 1/2. The foy (i ^ j) repre- 
sent cross-reactivities, which are generally beneficial for 
the immune system because they offer additional killing. 
However, the side effect of cross-reactivity is that it in- 
duces competition when the T cells expand in concentra- 
tion during the activation phase of the immune response. 

We first hypothesize an ideal immune system that has 
the following properties: (i) There is no competition 
among the different T cells. This allows for the free de- 
velopment of effector T cells whose total concentration 
will reach its maximum capacity, x = X)i=o x i = ^- For 
a realistic system, we will find x < 1 due to competi- 
tion, (ii) There is no bias toward any viral quasispecies 
in generating its corresponding effector cells. Therefore 
one has Xi/x = Vi/v. Since x = 1, one has Xi — Vi/v. 

The real immune system and the competition induced 
in the response to HIV-1 are modeled as follows. Com- 
pared to cytotoxic T cell lymphocytes (CTLs), antibodies 
are much less effective in controlling HIV-1. Therefore, 
we focus on the CTL responses, and a virus is accordingly 
identified by its epitope that stimulates CTLs. The epi- 
tope is N = 8-11 amino acids long, and in our model we 
assume N = 9. The initial strain of virus is denoted by 
Vo- The assumption of an unique initial strain is based 
on the bottleneck effect and the reversion effect due 
to lack of CTL pressure at infection ^| • We assume that 
at each of the N sites on the epitope there are A possible 



amino acid substitutions. The quasispecies space size is 
thus I = (A + 1) N . We use A = 2 because each site can 
accept few substitutions to keep the virus viable. We find 
that the simulation results are not sensitive to A. The 
subscript of Vi is in the range < i < (1 + A) — 1 
and can be written in form of a A + 1-radix number 
i = (ajy, cln— l) ' ' ' j with < a n < A. The Hamming 
distance between epitopes is = Yl n =id ~ ^a* a j )' 
where 5i m is the Kronecker delta. The model is 

dv ■ 

-T£ = TiVi - mNAvi + mTi - ci/^x)^ (1) 
-7^ = A - x^j - c 3 gi(x)Xi, (2) 

for i = 0, 1, • ■ ■ / — 1, with the initial condition vo(0) ^ 
0, Vi^o(0) — 0, and Xi(0) = 0. Previous models have 
assumed linear suppression (/; = const) and not tracked 
both viral and T cell quasispecies evolution 

In Eq. (Q, r,j — r exp(— adio) is the replication 
rate of Vi. The wild type Vo replicates the fastest, 
and mutation away demands some fitness cost. The 
term m is the rate of mutation of the amino acids 
at any one of the N sites into any one of the other 
A amino acids (per amino acid per day). The term 
mNAvi thus represents mutation away from qua- 
sispecies i to its NA neighbors (efflux). The term 
mTi represents the reverse process (influx), where Ti = 

X/j=0 ( W (i,nv-i,---,oi) + v (o,n ,j,—,ai) + • • • + v (aN,aN-i,---,j)) 

Nvi is the sum over all viral species a unit Hamming 
distance away from i. The term /i(x) = j/i(x)/[c2+yi(x)] 
represents the lysis ratio of targeted cells harboring 
Vi, where ^(x) = V' x,l-,, = Y^jZo ^j^Pi.- 
The derivation of /i(x) follows as in 13] by treat- 
ing the competitive binding process as a Langmuir 
adsorption isotherm of the CTLs onto targeted cells 
that are infected with Vi. The concentration and 
affinity here are dimensionless, while in [lflj their 
physical values {xj} and {kji} are used. That is, 
/((x) = ^ {.XV }{%}/[! + V ,{.r. ! {/,-.,} . We let 
{zjmax and {fcjmax denote the maximum values of 
concentration and affinity, respectively. We, thus, obtain 
c 2 = l/({a ; }max{^}max)- Typical concentrations are on 
the order of 10~ 6 mol/1, typical binding constants are 
on the order of 10 6 1/mol, and so we expect c 2 to be on 
the order of unity. Indeed, to fit clinical data, we will 
find C2 = 0.8. The coefficient c\ is the rate of clearance 
of virus released from targeted cells that are being lyscd, 
and to fit clinical data, we will find c\ = 4.0 day -1 , close 
to a previous estimation of 4.3 day -1 [l4j| . 

The first part of the right hand side of Eq. J2J repre- 
sents a delay process with time constant 1/A, since the 
immune system always lags in response to viral changes. 
The term Vi/v, which equals Xi in an ideal immune sys- 
tem, now serves as the moving target for Xi to catch. 
Moreover, the tracking of the virus by the immune sys- 
tem is disturbed by the competition c^gi(x)xi, where 
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<7i(x) = J2j^i x j^ji — Vii x ) ~ x i represents competition 
from the other T cells. This term has been used in study- 
ing the competition within an ecosystem . The T cells 
that develop into JQ, must be stimulated by V^. However, 
since Xj (J ^ i) can occupy Vi, a delay in the maturation 
of the T cell precursors into Xi can occur. 

The number of non-linear equations in the model is 
2(j4 + 1) N = 39366, which is difficult for analysis and 
computation. However, there is apparent symmetry in 
these equations: u, = Vj and Xi — Xj as long as dio = djo- 
That is, all the viruses that have an identical distance to 
the wild type have the same dynamics and can be put into 
the same group. We let v n (x n ) denote the number of 
elements in the viral (CTL) group n, one has v n — Q^Vi 
and x n — Q'^Xi for any i satisfying d^ = n, where — 
C%A n and C% = N\/[n\(N-n)\}. By substituting v n and 
x n into Eqs. and noting that = (n/Q^ 1 )v n -i + 

[n{A - 1)/Q%]v n + [(N - n)A/Q n N ' 

dv n 



" r ^Vn+it one obtains 



dt 



~dt~ 



mA (N — n + 1) v n —i + m(n + 1) v n 

+ (jo exp (-an) - c x f n (x)) v n 
—mA (N — n + n/A) v n 



+ i 



(3) 



v 



c 3 y n (x) x n + (4) 



for n = 0, 1, • • • ,N, where /„(x) = j/„(x)/[c 2 + 2m(x)]. 
To obtain y n (x) from yj(x), the following two steps are 
performed. The first step is to rearrange the equation for 
Hi by summing up all the xj that have the same d 



namely to obtain y„(x) = Yln'=0 1 n ' exp(— /3n'). To 



be specific, let us consider 
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with doi — n and a specific dji — n' , j must sat- 
isfy n — n' < djo < n + n' . To obtain a j with 
dji = n' and djo = fc, some 0's or l's in this expres- 
sion for i need to be changed, and the number of 
ways to realize this is Y%=k- n Q l N-n C n P t-h> where 
Q l N-n = C l N ^ n A l corresponds to changing I of the 0's 
into non-0; C^ 1 corresponds to flipping l\ of the l's 
into 0; P^_i ± = ^'^^(A — l)' 2 corresponds to changing 
I2 of the remaining (n — li) l's into digits other than 
and 1. The restrictions are I + h + I2 = n' and 
n + I — l\ = k, from which one solves l\ = n + I — k 
and I2 = n' — n + k — 21. Therefore one obtains q n i = 



E n + n (Z. m k \ X^ n D l C'n+l-k pn -n+fc-2f 

k=n—n' \ x k/ ^ N I Z^l=k—n "*JV— n r> r k-l 

The second step is to sum up all terms that have the 
same Xk- That is, to obtain the matrix that is 
defined by y = <f>x. By substituting the expression 
of q n i into the expression for y n one obtains = 

En'=oleM-(3n')/Q k N ] Etk-n Q l N- n C: +l - k P^7 n+k - 21 
The model parameters are determined as follows. The 
term mA is the mutation rate per amino acid site per day. 
The conventional HIV-1 mutation rate n w 2.2 x 10 -5 is 



per nucleotide base per replication 0|. Given that one 
replication cycle is about one day, one amino acid has 3 
nucleotide bases, and about 3/4 nucleotide mutations are 
nonsynonymous, one has m = 3/xx0.75/A w 2.475 x 10 -5 
per amino acid per day. The other parameter values are 
estimated by fitting the standard data given by 0. The 
viral load data is in terms of the plasma HIV-1 titer, 
which is two raised to the number of serial x 1/2 dilutions 
that result in undetectable virus. In the first vircmia 
phase, the mutation has not yet taken significant effect, 
Eqs. llll'll can be simplified as vq = vqVo — c\XqVq/(c2 + 
xq) and ±0 — A(l — xq). Since the CTL response peaks 
(i.e., xo a ppr oaching 1) at 9 to 12 weeks after the initial 
infection [ljj, we calculate from x — A(l — x Q ) that A 
is 0.033 day -1 . By noting that the first viremia peaks 
roughly at week 6 with plasma virus titer reaching 630 
and declines by week 12, we find «o(0) = 3.7 x 10~ 7 , ro = 
1.88 day" 1 , a = 4.0, and c 2 = 0.8. The fitting of the 
latency and the final viremia phases requires a — 0.41, 
/3 = 0.58, and C3 = 9.56 x 10~ 2 . As discussed above, only 
a and C3 lack comparison to literature values. 

Eqs. (jniEJl are still difficult for analysis because v shows 
up in the denominator of Eq. By summing Eq. Q 
over n, one obtains 

dv N 

~d7 = ^2 ( 7o ex P(~ an ) - Ci/„(x)J v n . (5) 

By defining p n — v n /v and plugging dv n /dt = vdp n /dt + 
p n dv/dt into Eq. |J5J, one obtains 



dpn 
dt 



N 

-Pn P" 
n'=0 



dX n 

~~dt 



mA(N - n+ l)p„_i + 

70 exp(-an) - ci/„(x)^ p n 
-mA(N — n + n/A)p n + rn(n + l)p„ + i 

„/ (j a cxp(-an') - ci/„'(x)) (6) 

A(p„ - x n ) - c 3 y„(x)i„ + -p^-Xn, (7) 

which are regular nonlinear ODEs. By integrating the 
above equations we find that p n (t) and x n (t) approach 
the fixed point x n ) as early as day 600. The eigenval- 
ues of the Jacobian matrix evaluated at the fixed point 
are all negative, indicating that the fixed point is sta- 
ble. That is, an equilibrium has been reached between 
the HIV-1 quasispecies and the immune system T cell 
quasispecies. This result holds even as the viral load it- 
self increases exponentially. All quasispecies escape at 
the same rate, with their relative proportions fixed. To 
calculate the rate of escape, Eq. (0 is converted into 
v(t) — e(t)v(t), where 



e(t) 



N 

^Pn ( 7 exp(-an) 

n=0 



Cl J2k ®nkXk 
C 2 + Dfc ®nkX k , 



(8) 
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FIG. 2: The escape rate e increases as the competition C3 
increases. 



Because the viral quasispecies converges, e(t) will even- 
tually converge to a constant e. Replacing (p n ,x n ) by 
the fixed value (p n ,x n ), one obtains s = 0.0066 day -1 . 
Because the equilibrium begins early, we conclude that 
the majority of latent phase, which was previously mys- 
terious, follows the same simple dynamics as the final 
viremia phase. This result explains the observation 0] 
that the lifespan of HIV patients is inversely correlated 
with the minimum viral load t> m i n that happens at about 
day 150 and is a direct result of the primary immune re- 
sponse. Because the onset time to the final phase scales 
as t oc — lnw m i n /e, the smallness of e significantly ampli- 
fies the life span dependence on u m i n . 

By incorporating immune competition, we now have a 
physical theory based upon an analytical model. Upon 
infection, the virus proliferates with a large rate around 
1.88 day -1 , because the immune response has not been 
fully established. At about week 6, the immune sys- 
tem brings the virus down to a trace level. However, 
the ever increasing viral genetic diversity induces com- 
petition among CTLs, which degrades the quality of im- 
mune control. Both a viral and a T cell quasispecies are 
established. The concentrations of different CTLs be- 
come fixed, as well as the percentages of different viral 
quasispecies. However, the virus escapes exponentially. 
Interestingly, it is this underlying pair of quasispecies 
that drives the seemingly nonequilibrium process. While 
small, C3 has become the Achilles' heel of the immune 
system. Competition among T cells may exert a small 
influence for many diseases, but is fatal for HIV-1. 

Our analysis predicts that competition among CTLs 
has a pivotal role on viral dynamics. The function 
e = e(c-z) is shown in Fig. One sees that as C3, or the 
level of competition, increases, the escape rate e increases 
significantly. For the ideal immune system, C3 = and 



e = —0.0065, suggesting that the ideal immune system 
can control HIV-1 infection. Even for the real immune 
system, the smallness of \e\ implies that HIV-1 and the 
immune system are well-matched competitors. When all 
antigens are simultaneously present, competition among 
T cells for activation and division occurs, and C3 is large. 
By vaccinating to different lymph nodes with the differ- 
ent HIV antigens likely to occur, competition between T 
cells is reduced, and T cells that recognize each of the 
different antigens may be produced at high concentra- 
tion 0], and so C3 w 0. In this case the virus never rises 
again after week 12. The regimen should be designed 
such that its components are given according to the fixed 
viral quasispecies. According to our simulation, the final 
viral percentages other than p and p^ are negligible. As 
for sensitivity to the parameters, if a or (3 are increased, 
in which case the data of Fig. ^ are n °t fit so well, pe 
also becomes non- negligible. We thus propose a vaccine 
regimen that includes not only the wild type but also 
the ps quasispecies. This vaccine must be given before 
the eventual evolution of the quasispecies viral strains to 
avoid immunodominance. If the vaccine were given after 
evolution of the p§ species, then competition would be 
unavoidable due to prior exposure of the immune system 
to the multiple viral strains. 
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